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Abstract 



We study the exclusive semileptonic B meson decays B~ D*^ i~ v 
and B^ D*^ v using data collected with the CLEO II detector at 
the Cornell Electron-positron Storage Ring (CESR). We present measure- 
ments of the branching fractions, B{B'^ £~ i>) = (0.5//oo)[4.49 ± 
0.32(stat.) ± 0.39(sys.)]% and B{B- D*^ £- v) = (0.5//+_)[5-13 ± 

0.54(stat.) zb 0.64(sys.)]%, where /oo and /_| are the neutral and charged 

B meson production fractions at the T(45) resonance, respectively. As- 
suming isospin invariance and taking the ratio of charged to neutral B me- 
son lifetimes measured at higher energy machines, we determine the ratio 
/+_//oo = 1.04 ± 0.13(stat.) ± 0.12(sys.) ± O.lO(lifetime); further assum- 
ing /_! + /oo = 1 we also determine the partial width r(5 — > D* i v) = 
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[29.9 ± 1.9(stat.) ± 2.7(sys.) ± 2.0(lifctimc)] ns-^ (independent of /+-//oo)- 
From this partial width we calculate B D* i v branching fractions that 

do not depend on /_| //oo nor the individual B lifetimes, but only on the 

charged to neutral B lifetime ratio. The product of the CKM matrix el- 
ement \Vch\ times the normalization of the decay form factor at the point 
of no recoil of the D* meson, T{y = 1), is determined from a linear fit to 
the combined differential decay rate of the exclusive B D* i v decays: 
\Vcb\^{^) = 0.0351 ± 0.0019(stat.) ± 0.0018(sys.) ± 0.0008(lifetime) . Using 
theoretical calculations of the form factor normalization we extract a value 
for iKbl- 
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I. INTRODUCTION 



In the framework of the Standard Model of weak interactions the elements of the 
3x3 Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix must be determined ex- 
perimentally. The element \Vcb\ is determined from studies of the semileptonic decays of 
B mesons. Measurements of \Vcb\ from the inclusive semileptonic rate 0,^3 from 
exclusive rates ll5|J^.|^,3flCT are systematically limited by model dependence in the the- 
oretical prediction of the decay rate. The recent development of Heavy Quark Effective 
Theory (HQET) ||lT| yields an expression for the B ^ D* i u [jl2| decay rate in terms 



of a single unknown form factor |TB| which, at the point of no recoil of the D* meson, is 



absolutely normalized up to corrections of order l/wig (where rriQ is the 6 or c quark 
mass). It is currently believed that these corrections can be calculated with less than 5% 
uncertainty |]15|,|16[, which would permit a precise determination of \Vcb\ from the study of 
B ^ D* i 1? a.s a. function of the recoil of the D* meson. The decay mode B —>■ D* i u 
is also preferred over other exclusive channels because D* meson decays have a very clean 
experimental signature. 

Throughout this paper the square of the four momentum transfer in B ^ D* i u decays 
is denoted by = M^^ (where Mm is the mass of the virtual W). The kinematic variable 
of HQET, which is a measure of the recoil of the D* meson, is given by 

where v and v' are the four velocities of the B and D* mesons and and m^). are their 
respective masses. 

We report on new measurements of the branching fractions and differential decay rate 
for the decays, 5° ^ D*+ t v and 5" ^ z/ [O. The CKM matrix element \Va>\ 



is extracted from fits to the differential decay rates as well as the integrated rate. Using 
isospin invariance to equate the partial widths of B^ D*^ l~ v and B^ —> D*^ i~ u, and 
recent B meson lifetime measurements, we also measure the T{AS) branching fractions, 

/oo = BiT{AS) ^ B^B') (2) 
/+„ = B{T{AS) B-B+). (3) 

All measurements of exclusive B branching fractions at the T{4S) resonance currently as- 
sume equal production of charged and neutral B mesons, so the measurement of these 
production fractions affects both hadronic and semileptonic B branching fractions. A model 
dependent result for the inclusive branching fraction B{B D* X i u), where X is any 
possible hadronic state (X 7^ 0), is also presented. 

The paper is structured in the following way: in Section II the technique for obtaining 
yields of B^ — > D*~^ i~ v and B~ — >■ D*^ l~ v events is discussed in general. Details of the 
candidate selection are then given in Section HI. Backgrounds in the B D* I v samples 
and how their magnitudes are estimated are discussed in Section IV. Systematic studies 
of the reconstruction efficiencies are described in section V. Branching fraction results are 
presented in Section VI, followed by measurements of |Vc6| in Section VII, and conclusions 
in Section VIII. 



4 



II. METHOD 



The — i> D*^ i~ V and B~ —>■ D*^ £~ u efficiency-corrected yields, denoted by A^^o and 
respectively, depend on the total number of T(45') decays in the data sample, iVx(4S), 
and a product of branching fractions, 

iVo = 4iVx(45)/ooi3(5° ^ D*+ r iy)BD*+BDO (4) 
N_ = ANri4S)U-l3iB- ^ D*"^ r z>)i3^*oS^o, (5) 



where Bd*+ = B{D*+ D°-k+), Bd*o = B{D*^ D^n^) and Bdo = B{D^ R- 



7r" 



These three branching fractions have been measured by CLEO II [|T8|,|T9[ with D* samples 
that are statistically independent of the sample considered here. The factors of 4 enter 
because each T(45') decay produces a BB meson pair and because we are combining the e 
and /i lepton species. 

We search for the decays B ^ D* i i?hj combining reconstructed D* mesons with "right 
sign" lepton candidates in the same event. By "right sign" we mean that a D*+ (D*^) must 
be paired with a negative lepton, while a D*~ requires a positive lepton. D* mesons 

are reconstructed using the decay chains D*^ D^7t~^, D*^ D^tt^ and K~7r~^ 
211 . This technique uses our knowledge of the B meson momentum, | p_B | . The energy of B 
mesons produced in symmetric e~^e~ annihilations, Eb, must equal the beam energy, which 
is precisely known from machine optics; hence \pb\ can be determined from Eb and the 



known B mass |21]. A kinematic constraint is obtained by writing the invariant mass of the 



emitted neutrino as 

pI = (Pb - Pd* - Pif (6) 

where p stands for the 4- vector of the particle in subscript. Expanding this equation results 
in 

pI = {Eb - EdhY - \pb\ ^ - IpdhI ^ + 2\pB\\PD'e\cosQ (7) 

where {Eb, Pb) is the B meson 4-momentum, {Eo^e, po'i) is the sum of the D* and lepton 
4-momenta, and G is the angle between the 3-momenta pD*e and p^. The ffist three terms 
on the right hand side constitute what is traditionally referred to as missing mass squared, 
symbolized MM^. The factor multiplying cosG will be denoted by C for cosine multiplier. 

MM^ = {Eb - ED*ef - \pb\ ' - \pD*e\ \ (8) 
C = 2|pb||p^*,|. (9) 

Ej:,*^ and PD*t are determined from the measured momenta of the lepton and D* candidates. 
Since we know the magnitude of p^, but not its direction, cosG is the only unknown. 

For each D* and lepton combination in the same event the pair of variables C and 
MM"^ is calculated. For correctly reconstructed B — »• D* I v decays (with perfect detector 
resolution), the values of C and MhP must lie within the kinematic boundary determined 
by Eq. with = 0. This boundary is shown in Fig. |I] for the lepton momentum 
range 1.4 < |p^| < 2.4 GeV. This lepton momentum range and kinematic boundary define 



the Signal Region ||22|. To arrive at the number of true B D* i v decays in the data 
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sample, we count the number of candidates observed in the Signal Region, subtract the 
expected number of background candidates which happen to fall inside the Signal Region, 
and divide by the Monte Carlo efficiency for signal events. Because the Signal Region spans 
a significant area of phase space, we cannot assume that backgrounds vary slowly in and near 
it, and a reliable estimate of the background inside the Signal Region cannot be obtained 
by interpolating from the number of candidates observed outside. Instead, we categorize all 
sources of background and estimate the total contribution from each source in the Signal 
Region by using data. The method is insensitive to the detailed C vs. MM"^ distribution of 
the signals and of all backgrounds except B —>■ D* X i u. 

B D* X i u decay is a significant source of background in this analysis, principally 
because this is the background physical process most similar to B ^ D* i v. It includes 
resonant B D** i v decays followed by D** D*X, as well as non resonant decays. 
This background is estimated from data events in a different region of C vs. MM"^ with 
0.8 < Ip^I < 1.4 GeV. This region in Ip^l, C and MM"^ is called the Correlated Background 
Region. It is described in more detail in Section IVB, along with the method for subtracting 
B D* X i u background. The method is sensitive to the detailed C vs. MM"^ distribution 
oi B —>■ D* X i u background, for which we use theoretical model predictions (and to which 
we assign conservative errors). Since this method must provide an estimate of the number of 
entries in the Signal Region due to B —> D* X i u decay, it also provides a model dependent 
measurement of B{B ^ D* X i u). 



III. EVENT RECONSTRUCTION 

The data used in this analysis were produced in symmetric electron-positron collisions at 
the Cornell Electron-positron Storage Ring (CESR) and recorded with the CLEO II detector. 
The signal comes from an integrated luminosity of 1.55 fb~^ collected at the T(4S') center- 



of-mass energy pSf. An additional 0.69 fb~ of data collected below the BB production 
threshold are used for continuum background determination. 

The most crucial components of the CLEO II detector in this analysis are the tracking 
system, the Csl electromagnetic calorimeter, and the muon identification system. A detailed 
description of the CLEO II detector is given elsewhere [0]. The tracking system comprises a 



set of drift and straw tube chambers in a 1.5 Tesla magnetic field that measure the momenta 
of stable charged particles over approximately 92% of 4tt with a transverse momentum 
resolution of [Spt/ptY = (0.0015pj)^ + (0.005)^, where pt is measured in GeV. Photons 
are detected in a Csl electromagnetic calorimeter with an angular acceptance of 95% of 
An. We restrict the fiducial volume for photons to the barrel portion of the calorimeter, 
I cos^-yl < 0.8, where 6.y is the angle a photon makes with the beam line (polar angle). The 
calorimeter energy resolution is AE/E{%) = 0.35/^^/^ + 1.9-O.lE, where E is in GeV p5|, 
which corresponds to 4% at 0.1 GeV. Both electron and muon candidates must lie within 
the polar angular region | cos^^l < 0.71. In this analysis, electrons with momenta above 0.8 
GeV are identified by their electromagnetic interactions in the calorimeter, their energy loss 
in the drift chamber gas and their time of flight in the detector. The electron identification 
efficiency within the fiducial volume is over 94%, while a hadron in the momentum range 0.8 
to 2.4 GeV has on average a (0.3 ± 0.1)% probability of being misidentified as an electron. 
Muons are identified by their ability to penetrate at least 5 nuclear absorption lengths in 
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iron, which puts a lower hmit of 1.4 GeV on the muon momentum acceptance. Muons 
within the acceptance are identified with 93% efficiency, while hadrons have on average a 
(1.4 ± 0.2)% probability of being misidentified as a muon. 



For this analysis we select hadronic events p6[ that have at least one track identified as 
a lepton with momentum 0.8 < |p^| < 2.4 GeV. The ratio of the second to the zeroth Fox- 
Wolfram moments of the event is required to be less than 0.4 to suppress background 
from continuum events. For each lepton in these events we search for candidates in the 
decay mode — > K~7r~^ using charge correlation with the lepton to make unambiguous 
mass assignments (the lepton and kaon charges must be the same). 

We combine candidates with pion candidates to fully reconstruct D* mesons in the 
modes D*^ D^n^ and D*^ — ^ D^tt^. We call these pion candidates slow pions because 
their momentum is restricted to be less than 225 MeV in the laboratory frame. Charged 
slow pions are accepted if they lie in the polar angle region | cos 6*7^+1 < 0.71 and have 
momentum above 65 MeV. Candidate vr^'s are constructed from pairs of showers in the 
electromagnetic calorimeter which do not match the projection of any drift chamber track 
and have an invariant mass within 3a of the measured 7r° mass (cr = 5 to 8 MeV, depending 
on shower energies and polar angles). Showers used in 7r° candidates must be in the polar 
angle region | cos^^| < 0.8 and have an energy above 30 MeV. The 7r° momentum vector is 
reconstructed by constraining the shower position and energy measurements to produce the 



known 7r° mass. The momentum of D* candidates must satisfy Ipd*]/ y-E's — "^d* < 0.5 to 
be consistent with B decay. 

The raw yield of events with a D* meson and a lepton is obtained by fitting the 
mass peak after cutting on the D* — candidate mass difference, 5m = Mk-ktt^ — Mxn, and 
subtracting the scaled result of a mass fft to a 6m sideband {Mxn and Mxtt-ks denote the 
invariant masses of the and D* candidates, respectively). This sideband subtraction is 
described in detail in the next section. 



IV. BACKGROUNDS 

The number of D* and lepton pairs observed in a given region of the C vs. MM"^ plane 
is the sum of the signal and background sources in that region. The background sources 
fall into ffve categories, and we have evaluated the contribution from each of these sources 
using the data. In order of importance the ffve background categories are: combinatoric, 
correlated, uncorrelated, continuum, and fake lepton background. 



A. Combinatoric Background 

The dominant background in this analysis is the combinatoric background in D*^ and 
D*^ reconstruction. One class of background arises from combinations of random K~tt^ 
pairs with any slow second pion, tt^. This class of background does not peak in Mk-k, but a 
subclass of these events in which the vr^ and either the K~ or 7r+ are from a true D* 
decay does peak in 5m- Combinations of K~tt^ pairs from a correctly reconstructed 
meson with an unrelated pion candidate forms a second class of combinatoric background 
which peaks in M^jr, but not in 5m- The Mk-k distributions in the 5m signal and sideband 
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regions are fit to a D'^ peak plus a Chebyshev polynomial, which removes the first class of 
background in both regions. To remove the second class and obtain a final signal yield, 
the D° yield in the 6^. sideband is scaled and subtracted from the yield in the 6m signal 
region. The scale factor is determined from the background function in a fit to the 



distribution 28 



The distributions of 5m for D*~^ and D*'^ are shown in Fig. ^ for events which lie in 
the C vs. MM^ Signal Re gion and have Mk-k within 100 MeV of the measured D mass. 
The 6m signal region for D*~^ is 8 MeV wide and is centered on the measured 6m mean of 
145.44 MeV while the sideband region is 150 < 6m < 165 MeV. The 6m signal region 
for D*^ is 6 MeV wide and is centered on the measured 6m mean of 142.12 MeV |3y], while 
the sideband region is 147 < 6m < 162 MeV. We have used a wider signal region for D*~^ 
because the 6m signal in this mode has non Gaussian tails from systematic effects in the 
reconstruction of very low momentum charged tracks. These effects are understood and 
reproduced in the Monte Carlo simulation. 

The Mktt distributions for events in the 6m signal and sideband regions are shown in 
Fig. ^. These are the four distributions that were fit to obtain the yields in the C vs. MM"^ 
Signal Region. The yields are listed separately for the 6m signal and sideband regions in 
Tables | and ^. 



B. Correlated lepton Background 

A B meson can decay to a final state with a D* and an e~ or /i~ through channels 
other than B —>■ D* i u, and these physical processes contribute a background that we call 
correlated. The principal source of correlated background is the decay B ^ D* X i i?, where 
e.g., X = vr or other unreconstructed particle(s) (X 7^ 0). To remove this background, we 
exploit two general features oi B ^ D* X i u decays. First, because we do not reconstruct 
X the MM"^ variable will tend to be shifted towards positive values. Second, because the 
D*X invariant mass is larger than the D* mass, the lepton spectrum will be softer than for 
B ^ D* i u decays. Fig. ^ shows the C vs. MM"^ distribution in the two lepton momentum 
ranges 0.8 < \pi\ < 1.4 GeV and 1.4 < |p£| < 2.4 GeV for Monte Carlo B D** i i> decays 
generated according to the model of ISGW |^ (we take these events to be representative, 
at the level of precision required here, of generic B ^ D* X I v decays). In the lower 
lepton momentum range there is a region of the C vs. MM"^ plane where the B D** I v 
efficiency is high, but there is almost no contribution from B —>■ D* i v decays. Therefore, 
the low lepton momentum range is used to measure the level oi B ^ D* X i u, which is 
then scaled using Monte Carlo to determine the B ^ D* X i u contribution to the high 



lepton momentum Signal Region [^. We do not determine the level of 5 ^ D* X i u 
from the events just outside the high momentum Signal Region, because the B D** i v 
efficiency is low in the high lepton momentum range and electron bremsstrahlung in the 
detector causes a small fraction of S — D* C. v decays to be reconstructed outside the signal 
boundary (Fig. |I|). 

The dashed lines in Fig. ^ define a region in the C vs. MM"^ plane for the lepton 
momentum range 0.8 < |p£| < 1.4 GeV called the Correlated Background Region, because 
it has high efficiency for the B D* X C v background processes, but very low efficiency for 
our signal mode. The M^tt distributions in the Correlated Background Region are shown 



8 



in Fig. ^. Let Ns and Nc denote the total yields in the high lepton momentum Signal 
Region and the low lepton momentum Correlated Background Region, respectively, after 
subtraction of all other backgrounds as discussed in the following sections. They are related 
to Ns{B — i> D* i u), the number of B ^ D* i v events in the high momentum Signal 
Region, and Nc{B D* X i u), the number of B ^ D* X i u events in the low momentum 
Correlated Background Region, via 



Ns = Ns{B ^ D* tu) + rcNc{B ^ D* X i u) 
Nc = rsNs{B ^ D* iu) + Nc{B ^ D* X i v) 



(10) 



where = 0.015 ± 0.003 is the ratio of the efficiency for B D* C. v events in the 
Correlated Background Region to the efficiency in the Signal Region, and Vc = 0.70 ± 0.35 
is the ratio of the B ^ D* X i u efficiency in the signal region to the efficiency in the 
Correlated Background Region. Even though the Correlated Background Region is outside 
the nominal kinematic boundary for signal, is non zero because of electron bremsstrahlung. 
The values for and Vc were determined from Monte Carlo simulation, and the errors reflect 
uncertainties due to model dependence. We estimated Vc using resonant B —>■ D** i v decays. 
These were generated according to the ISGW model pTI for the ^Pi, ^Pi, and ^P^ D** states 



in their predicted relative abundances We expect that the B D**{^Pi)lv decay is a 
reasonable approximation to non-resonant B —>■ {D*TT)iu decays because the ^Pi is expected 
to be a very wide resonance. We take into account our rough modelling of P — * D* X i u 
decays by assigning a conservative error of 50% to Tc. 

The observed numbers of events and the yields obtained by solving Eq. (|l^) with the 
quoted central values of and are summarized in Tables | and ||. The uncertainties in 
Tg and Tc are included in the systematic errors for the final B D* £ v and B ^ D* X i u 
yields. 

Sources other than B D* X i u contribute to the correlated background at much 
smaller levels. The decays B D*X, where X fragments to light hadrons that decay 
semileptonically or are misidentified as leptons, are discussed in the section on lepton fakes. 
The processes B D*t^u followed by l^uv , and B P*P|*j followed by 

p|*j^ Xl~u lead to final states that include a D* and a lepton with the same "right sign" 
charge correlation as signal. These sources have characteristics similar io B ^ D* X I 
but with an even softer lepton spectrum. Monte Carlo simulation predicts values for of 
0.01 ± 0.01 and 0.005 ± 0.005 for B D*t-v and B D*d\*J~ respectively. They 
can therefore be neglected in the Signal Region. We have estimated their contributions to 
the background region using available measurements or estimates of the various intermediate 
branching fractions p9| , |34| together with efficiencies determined in Monte Carlo simulations. 



The results are given in Tables | and ||. 



C. Uncorrelated D* lepton background 

"Uncorrelated" background encompasses events with a D* from the decay of the B and 
a lepton from the P in a PP event. However, in order to contribute to the background in 
this analysis the D* and the lepton must have the same charge correlation as signal decays. 
Therefore, the leptons in uncorrelated background are secondary P meson decay products 
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from the decay chain h c C (known as cascades) or ^ ^ iX (mixing). It is also 
possible to have leptons from the decay or misidentification of light hadrons {K or vr). 

The cross section for producing uncorrelated background events is measured using the 
data. Because the B and B mesons in an event decay independently, the uncorrelated 
background cross section is a product of the inclusive D* cross section from T{4S) decay, 
(T'{e~^e~ —>■ T{AS) D*), and the B ^ i~ inclusive branching fraction, B'{B The 
primes indicate that the quantities of interest are the detected ("raw") cross section and 
branching fraction, which include detector acceptance and reconstruction efficiencies as well 
as the underlying cross section and branching fraction. We measure these two quantities 
with the same event selection and cuts used to obtain B ^ D* i u yields. 

We find that the cross sections for producing and detecting a D* from T(4S') decay 
(integrated over momentum) are 3.1±0.15 pb and 3.2±0.33 pb for D*~^ and D*^ respectively. 

Measuring B'{B — ^ i~) is not as straightforward as counting leptons because of the 
need to preserve the correct D* and lepton charge correlation, i.e., we need specifically 
B'{B — >• rather than B'(T{4:S) — >• £^). We measure B'{B — >• i~) using like charge 
dilepton events, in which one of the leptons tags one daughter of the T{4S) as either a B 
or a 5 |35|. We find that B'{B £-) is (0.7 ± 0.05)% and (1.2 ± 0.05)% in the momentum 
range 1.4 < Ip^l < 2.4 and 0.8 < |p^| < 1.4 GeV, respectively. These detected branching 
fractions include many sources, such as misidentified hadrons from the "other" B, photon 
conversions, etc., but these all contribute to uncorrelated background and therefore belong 
in the quantity we are trying to measure. 

Multiplying a'(e+e~ T{AS) D*) by B'{B £') yields the total cross section 
for uncorrelated background. However, we wish to estimate the number of uncorrelated 
background events in a particular region of the C vs. MM"^ plane; therefore the distribution 
of uncorrelated background in these variables must be known. These variables depend only 
on the inclusive D* and lepton momentum distributions (which we have already measured) 
and the angle between the D* and the lepton, a. The distribution of cos a is fiat because 
the two B mesons are produced nearly at rest. We use these three known distributions to 
simulate the C vs. MM"^ distribution of uncorrelated background. From this simulation and 
the measured uncorrelated background cross section we estimate the numbers of uncorrelated 
background events in the Signal Region and Correlated Background Region (given in Tables 
|and[T|). 

D. Continuum and Lepton Fakes 

The level of continuum background is estimated using the 0.69 fb~^ of data recorded at 
energies slightly below the T(4S') resonance. These data are analyzed in the same manner 
as resonance data in order to estimate the continuum backgrounds listed in Tables | and ||. 

There is a small background of D*£ candidates where the lepton is actually a hadron that 
has been misidentified as a lepton. We call this the fake lepton background. Fake D*i pairs 
are predominantly uncorrelated, simply because there is more energy available to produce 
light hadrons in the decay of the B meson that did not produce the D*. These uncorrelated 
fake leptons are already included in the cross section for uncorrelated background. 

The correlated lepton fake contribution to the raw number of D*l pairs is estimated by 
performing a similar analysis of the data where D* candidates are paired with light hadron 
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candidates instead of leptons, and scaling the result by the known electron and muon fake 
rates. For this study a light hadron is defined as any detected charged particle that fails very 
loose electron and muon identification cuts. We use Monte Carlo simulation to determine 
the efficiency for correlated D* and hadron pairs to fall within the signal and background 
regions of the C vs. MM'^ plane. The correlated lepton fake contributions turn out to be 
so small that we can easily afford a large uncertainty due to the use of Monte Carlo. The 
results are given in Tables | and y. 



V. DETECTION EFFICIENCIES 

The efficiency for B —>■ D* i i? events to pass our event selection criteria is estimated from 
a Monte Carlo simulation. The generator produces events obeying all kinematic constraints 
due to angular momentum conservation and the V — A nature of the pseudoscalar to vector 
decay [^. These events are passed through a full detector simulation and then analyzed 
as described above. 

An uncertainty in the efficiency for B ^ D* i u exists because the dynamics of the 
decay depend on unknown form factors [|38|. Several phenomenological models of these 



form factors are available. For the results presented here, we have used the predictions of 
Neubert |3^;^^ to determine the central values of our efficiencies, and have compared these 



results to those obtained with the ISGW BSW [^], and KS models to estimate 



the model uncertainty (See Tables |V| and |V11| ). After comparing these models as well as 



varying the form factors within the Neubert model, it is estimated that the model dependence 
contributes a 3% systematic uncertainty in the detection efficiency. 
The efficiencies for B D* I v events to be counted as signal are 

e,{& D*+ r v) = [9.54 ± 0.23 ± 0.71]% (11) 
esiB- r u) = [7.18 ± 0.52 ± 0.53]% (12) 

where the first (systematic) error is correlated between the two efficiencies and the second 
is uncorrelated. These efficiencies do not include the D* and D'^ branching fractions. The 
errors include: Monte Carlo statistics, model dependence, an uncertainty reflecting changes 
in efficiency corrected yields due to variation of cuts, the uncertainty in the normalization 
of the 6m sidebands, the uncertainty in the lepton reconstruction and identification, a 2% 
per track uncertainty in the reconstruction of K and tt tracks with momenta above 250 
MeV, and the uncertainty in modelling the efficiency of the slow pions. The methods used 
to estimate the latter uncertainty are described below. The various contributions to the 
systematic uncertainty in the efficiency are listed in Table fT|. 

In the laboratory frame, the slow pion from the B D* i u, D* Dtt decay chain 
has a momentum of less than 225 MeV with the peak of the distribution at 100 MeV. The 
efficiency to reconstruct vr'^'s is fairly fiat in this momentum range but difficult to model 
because of the large background of low energy showers in the calorimeter. The efficiency for 
reconstructing vr+'s decreases sharply below 100 MeV, falling to at 50 MeV. It is crucial 
to this analysis that these efficiencies be correctly reproduced by the Monte Carlo and we 



have performed several detailed studies to evaluate the reliability of our simulation [43 
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A. vr"'" Efficiency 

The shape of the efficiency curve for charged pions of momentum less than 225 MeV can 
be measured using inclusive D*^ decays in the data. In the strong decay D*^ —>■ D^7r~^, 
dN / d cos if D'-K niust be symmetric about cosipD*.,, = 0, where (pD*n is the angle between the 
slow vr'^ momentum in the D*~^ rest frame and the D*~^ direction in the lab frame. However, 
the observed dN/d cos ipD*n distribution in the data may be asymmetric because the slow 
TT"*" efficiency varies with pion momentum, which is highly correlated with cos ^. 

The charged pion efficiency between and 225 MeV was measured in data by simultane- 
ously fitting the observed dN / d cos (pn*.,^ distributions in 8 bins of D*^ momentum between 
and 5 GeV. It was found that the Monte Carlo reproduces the data efficiency shape over 
the pion momentum range to 225 MeV (Fig. |^). This study of the cosyj^i*^ distribution 
tests the simulation of the shape of the efficiency curve, but not its absolute normalization. 
We test the normalization by comparing the yield of fully reconstructed — >■ K~t:^t:^ 
to the yield of partially reconstructed —>■ K^7i^{7i^), where the 7r+ is not used. The 
momentum range for these pions is 0.2 to 1.0 GeV. The Monte Carlo efficiency agrees well 
with the absolute tt"'" efficiency measured in the data, within the statistical precision of 2%. 
Together, these studies check at the 5% level the simulation of the slow 7r+ efficiency. 



B. tt^/tt^ Efficiency 

The accuracy of the simulation of the slow 7r° efficiency is difficult to test. The method 
outlined above cannot be used for the decay D*^ — > D^tt^ because of the large combina- 
toric background in the vr*^ reconstruction. Instead, we compare the ratio of efficiencies, 
e{n'^)/e{TT~^), in the pion momentum range to 225 MeV between data and Monte Carlo. 
We then combine the results of this study with the uncertainty in the charged pion absolute 
efficiency described above in order to obtain the uncertainty in the tc^ absolute efficiency. 

We determine the ratio of efficiencies below 225 MeV, ci = e(7r'^)/e(7r+), using a combina- 
tion of two methods. The ffist part uses the decay 1] tt+tt^tt^ as a source of both neutral 
and charged soft pions. Candidate r/'s are selected with either a tt*^ in the low (0 to 225 MeV) 
momentum range and both charged pions above 250 MeV or a tt^ in the low momentum 
range and the other charged and neutral pions with momentum above 250 MeV. The ratio 
of the number of events in these two bins depends on the efficiency to find the pions times 
a kinematic factor that depends on the rj Dalitz plot and the t] momentum spectrum. The 
kinematic factor can be interpreted as the average probability for the three pions from the 
f] decay to populate the desired momentum ranges. The r] Dalitz plot has been precisely 



studied in previous experiments and the rj momentum spectrum in the Monte Carlo 
was tuned to the data using the 77 — > 77 decay mode. The ratio of yields is proportional 
to the ratio Ci/c2 where Ci is defined above and C2 is the ratio of neutral to charged pion 
efficiencies above 250 MeV. The ratio Ci/c2 is determined in data and compared with the 
same ratio calculated from Monte Carlo. 

A second study using Kg^s then determines the product C1C2. The Kg study takes 
advantage of the fact that the ratio of branching fractions B{Kg -k^-k^) / B{K^ tt+tt^) 
is well known. Therefore, the ratio of neutral to charged final state yields is the product of 
a known ratio of branching fractions, an efficiency ratio and an acceptance factor which is 
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determined from Monte Carlo. For kaons of momentum 250 to 500 MeV, one daughter pion 
is in the range 60 to 225 MeV and the other daughter is above 225 MeV so that the ratio 
of yields N{K'^ 7r°7i°)/N{K^ Tr+Tr") in the i^^ momentum range 250 to 500 MeV is 
proportional to the product C1C2. 

When the results of the t] and Kg studies are combined we find that ci (the ratio of 
neutral to charged pion efficiency for momentum 60 to 225 MeV), as determined from the 
data, agrees with the predictions of the Monte Carlo to within the statistical precision of 
7%. We therefore assign a systematic error on the efficiency ratio of neutral to charged slow 
pions of 7%. Combining this with our study of the absolute slow charged pion efficiency, 
we assign an error of 8.6% to the absolute determination of the neutral pion efficiency 
in the momentum range below 225 MeV. The Monte Carlo prediction of the efficiency to 
reconstruct neutral pions is shown in Fig. |^. 



VI. BRANCHING FRACTION RESULTS 

In this analysis we measure two independent quantities: A^o and of Eqs. (^) and (^. 
These are in turn given by 

where Ns{B^ D*^ i~ v) and Ns{B~ D*^ i~ u) are the background subtracted yields 
for the given signal modes (Tables | and|I|), and 6^(5° D*^ £~ v) and es{B~ D*^ i~ v) 
are the efficiencies for the given signal modes (Section [V]). Dividing A^o and by four times 
the number of T(4S') decays in the data sample yields the product branching fractions given 
in Table |V[ The number of T(4S') events was determined by studying the hadronic cross 
section at energies both on the T(4S') resonance and slightly below the resonance in the 
continuum. For our data sample, A^t{45) = (1-65 ± 0.01) x 10^ events. 

The top two product branching fractions of Table |V| depend on the unknowns to be 
determined, E{B^ D*^ i~ u) and B{B^ D*^ i~ i?), as well as the two unmeasured 
quantities /+_ and /qq. In all, there are four unknowns and only two measurements. The 
assumption that the T(4S') decays only to BB (/+_ + /oo = 1) reduces the number of 
unknowns to three. This assumption is supported by the observation that the T(4S') width 
is three orders of magnitude larger than lower T states which are not massive enough to 
decay to a pair of B mesons There are two ways to further reduce the number of 



unknowns. The first method is to make the traditional assumption that = /oo, which is 

uncertain at the 5 — 10% level Using this assumption we present the two independent 



measurements: B{B^ — >■ D*~^ i u) and B{B —>■ D*^ i u). Alternatively, we can assume 
that B{B- D*^ £- p)/B{B^ D*+ u) = tb-/tbo. This is equivalent to assuming 
that the partial widths for the exclusive semi-leptonic decays of the charged and neutral 
B's are equal: T{B- i~ v) = T{B° D*+ t' v). With this assumption and 

a measurement of Tb-Itj^q from other experiments our measurement of Nq and 
can be used to determine V{B D* i u) and /+_//oo. We will first discuss the B —>■ 
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D* i u branching fractions obtained using the two different assumptions, and then the 
T{AS) branching fractions obtained with the second assumption. 



A. B ^ D* i u Branching Fractions 



To extract the B meson semileptonic branching fractions from Table |I^ we require values 
for the rest of the branching fractions in each product. For the D* branching fractions we 
use the CLEO II results 



B{D*+ 



[68.1 ±1.0 ±1.3]% 
[63.6 ± 2.3 ± 3.3]% 



For the D branching fraction we use the CLEO II result [19]: 



[3.91 ±0.08 ±0.07 ±0.16]% 



(15) 



(16) 



which is the value given in Ref. [jl9| without final state radiation. The systematic error on K 
and TT track reconstruction (third error) has been separated from other systematics (second 
error) because it is correlated with row 5 of Table |ITT[ This correlation is an important 
advantage of using this CLEO II measurement of B{D^ K^tt^), as the error due to Ktt 
track reconstruction will cancel in the determination of the B D* i v branching fractions. 
Thus the fractional systematic error on our results for BiB D* I v) will be reduced 



relative to the error on the product branching fractions listed in Table |IV. 

If we assume equal production of charged and neutral meson pairs at the T(4S'), /4 
/oo = 0.5, we obtain. 



^0.5(5° 
^0.5(5- 



[4.49 ±0.32 ±0.39]%, 
[5.13 ±0.54 ±0.64]%, 



(17) 
(18) 



where the first error is statistical, and the second (systematic) error includes the errors on 



the and D* branching fractions. On the other hand, by assuming Yilf' 
Y{B~ —>■ D*^ i~ u), which follows from isospin invariance, and /+_ ± /oo = 
(H) can be combined to give 



Eqs. 



T{B D* 



1 



No 



+ 



N_ 



4A^T(45)^D0 [tboB£)*+ Tb-Bd*o 

Using the average of LEP and CDF measurements for the lifetimes tb- and r^o 
V{B ^ D* iu) = [29.9 ± 1.9 ± 2.7 ± 2.0] ns~\ 



§ and 

(19) 

yields, 
(20) 



independent of /+_//oo. The second (systematic) error takes into account correlations be- 
tween the errors of Eqs. ( p!7D and (|18]). The third error is due to the B lifetime measurements 
and is determined with the conservative assumption that the errors on the charged and neu- 
tral lifetimes are fully correlated. 

The partial width can be converted to either the charged or the neutral B —>■ D* i u 
branching fraction by multiplying Eq. ([19|) times the appropriate lifetime, but note that 
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B{B- 
result is 

B{B- 



D 



V 



and BiB"^ 



are not measured independently this way. The 



D 



*0 /)- 



BoAB' 



D 



*+ p- 



BoAB- 



D 



*0 /)- 



(21) 



independent of /+„//oo. Note that this does not depend on the individual B lifetimes, but 
only on their ratio. Taking Tb- /t^o = 1.1 ± 0.11 from other experiments leads to 



B{B- 



D 



*o 



1.1 X B{B^ D 



*+ 



[5.03 ±0.32 ±0.45 ±0.24]%, (22) 



where the uncertainty in our knowledge of the lifetime ratio leads to the third systematic 
error. Use of a different central value for the lifetime ratio will change the right hand side 
of Eq. ( p2D according to Eq. (|2l|) . The factor 1.1 carries no error to indicate that we have 
assumed that the ratio of branching fractions exactly equals the lifetime ratio, which means 
that we can not independently determine the charged and neutral branching fractions. 

One can compare the above branching fractions to previous results from ARGUS and 
CLEO, which assume /+_ = /oo = 0.5, after correcting all results for the new CLEO 11 D* 
and branching fractions The results of this comparison are shown in Table |V[ Our 
results are in good agreement with previous measurements. 



B. T{4S) Branching Fractions 

By taking the ratio of Eq. over and assuming only that r(i?° D*^ v) = 
T{B- i- z/), we have 

= = 1.14 ± 0.14 ± 0.13. (23) 

where the first error is statistical and the second (systematic) error is dominated by the 
uncertainty on the slow pion efficiency ratio, e,ro/e7r+- Substituting the value for tb-It^o 



from Ref. 47 



^ = 1.04 ±0.13 ±0.12 ±0.10, (24) 
/oo 

where the third error is due to the input lifetime ratio. Finally, the assumption /+_ ±/oo = 1 
can be used to extract a value for /+_ or /oo, 

/+_ = 1 - /oo = 0.510 ± 0.052, (25) 

where the error is combined statistical and systematic, including the error in the B lifetime 
ratio. 
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C. B{B D* X I v) from Correlated Background Region Yields 

The product branching fractions for B{B D* X i i?) in Table |V| have been computed 
from the correlated background yields in the Correlated Background Region from Tables | 
and |l| divided by an efficiency estimated with Monte Carlo B D** i v events, as explained 
below. 

All possible resonant and non resonant sources of correlated background contribute to 
the data yields: 

N^{B ^D* X iv)=MY. B{B D** i v)B{D** D*X)ei 

i 

+ AfJ2 D*Xj i u)ej, (26) 

j 

where i ranges over all possible resonant states and j over all non resonant states, and A/" 
contains the total number of B decays in the data times the D* and branching fractions. 
The efficiencies for the individual channels, and e^, are potentially all different. Therefore, 
an accurate estimate of the the efficiency for the B ^ D* X £ u inclusive process, e^B — >■ 
D* X i 1^), requires knowledge of the relative abundances of all the resonant and non resonant 
exclusive modes. However, if the different exclusive modes all have similar efficiencies, one 
need only simulate one or a few of them in order to estimate ec{B —>■ D* X i u). We 
compared the Monte Carlo efficiencies for B D** i v where the D** was a l^Pi, \^P\ or 
1^P2 state generated according to the ISGW model [pl|] in the D*7r decay channel. Since 



the l^Pi resonance is very wide, we expect it to resemble non resonant decays where Xj is a 
pion. For these three exclusive channels we find the following efficiencies in the Correlated 
Background Region: (3.4 ± 0.2)%(liPi), (2.3 ± 0.2)%(l3Pi), and (2.9 ± 0.2)%(l3p2) for a 
in the final state, and (2.5 ± 0.2)%(liPi), (1.7 ± 0.2)%(l3Pi), and (2.3 ± 0.2)%(1='P2) 
for a D*^ in the final state. We take the average of these numbers as a rough estimate of the 
Correlated Background Region efficiency for any channel which leads to a in the final 
state, 

ec(P ^ X r v) = [2.9 ± 0.6]% (27) 
e^B P)*° X £- iy) = [2.2 ± 0.4]%, (28) 

where the errors were estimated from the spread in the efficiencies of the three individual 
modes considered, together with the uncertainties from Table |T|. 

Dividing rows 3 and 4 of Table |^ by the and D* branching fractions, we obtain, 

B{B D*+ X r u) = [0.6 ± 0.3 ± 0.1]% (29) 
B{B P)*° X £- u) = [0.6 ± 0.6 ± 0.1]% (30) 

These inclusive branching fractions do not depend on the unknown relative abundances 
of the different resonant and non resonant D*X states in semileptonic B decay, as long 
as they all have similar efficiencies in the Correlated Background Region. To obtain an 
upper limit for B{B D**iu), we assume that non resonant channels are all zero, and 
use the relative abundances of the first radial excitation D** states in the ISGW model to 
estimate B{D** D*n) = 77% H. Together with resuhs (|2|) and (^) this leads to 
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EiB{B ^ D**ii?) = (1.5 ± 0.8 ± 0.2)%, which does not include an error on the ISGW 
estimate of the relative D** abundances. Converting this to an upper limit yields, 

^ 8{B D**iij) < 2.8% (31) 

i 

at the 95% confidence level. This result is consistent with previous model-dependent deter- 
minations of the B — s> D** i V branching fraction [^j. 



VII. DETERMINATION OF \Vcb\ 



The decay amplitude oi B ^ D* i v for massless leptons |^ is commonly expressed 
in terms of the three meson form factors A\[(^\ ^2(0'^) and Vis^^ These form fac- 

tors characterize the transition between the two strongly bound states B and D*, and are 
therefore not calculable. In order to calculate the decay rate for B ^ D* I v and hence de- 
termine \ych\ from the measured rate it is necessary to know the normalization as well as the 

dependence of these form factors. Several phenomenological models have been devised in 
order to estimate the form factors, but these typically assume some dependence and use 
ad hoc wave function representations for the B and D* bound states. Hence it is difficult to 
estimate the accuracy of model predictions. In order to compare different models we define 
the form factor ratios, AijAx and V IA\. Table [V| lists the expressions for Ai, A^jAx^ and 
V I A\ for the different models examined. 

Recent work on Heavy Quark Effective Theory and its application to the decay B — 
D* I V has led to constraints on the form factors that permit a less model dependent 
determination of |\4f,|. In the limit of infinite h and c quark masses the meson form factors 
A\^ A2 and V are well defined functions of a single form factor C,{y), known as the Isgur-Wise 
function [|ll|], where y is the kinematic variable of HQET, and is related to by Eq. The 



Isgur-Wise function ^{y) contains the non perturbative QCD dynamics of the light degrees 
of freedom in the mesons, and it is therefore not calculable. However, at y = 1 (point of 
maximum g^) it is normalized to unity. 

For finite b and c quark masses the B ^ D* i u differential decay rate can be expressed 
in terms of a single unknown form factor J-'iy) which incorporates the three form factors Ai, 
A2, and V. The relationships between these and the Isgur-Wise function are given by HQET 
only to leading order, and the corrections necessary to account for the deviation from the 



heavy quark symmetry limit are absorbed into J^iy). Following Neubert |]39| , |50[| we write 



riv 



v' - 1 



2 



(32) 



where r = m£,*/mB. The function J-'{y) can be related to the Isgur-Wise function and 
correction terms that vanish in the infinite b and c mass limit, J^(l) = ?7^^(1) + O (^^' 



QCD ' ^ 



where rjA is a perturbatively calculable QCD radiative correction and Aqcd is the QCD scale 
parameter. A next to leading order calculation gives tja = 0.986±0.006 ^9|]. For B ^ D* i P 



decay, it has been shown that corrections of order l/rric and 1/mi, are identically zero [jT4 
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and it is currently estimated that the second order deviations from the symmetry hmit can be 



calculated with approximately 4% uncertainty ||15| , p!6| , |39[| . Therefore, a precise measurement 
of |V^f,|jF(l) will result in an accurate determination of |Vcf,| with no need for predictions of 
the shape of the decay form factors. Measurements of the differential decay rate dT / dy will 
also determine the shape of J-'{y), which provides information about the non perturbative 
QCD dynamics of the decay. 

In the next section we use model predictions for the Ai, A2 and V form factors to 
extract model dependent values of \Vcb\ from our measurement of B{B^ D*^ l~ u) and 
E{B^ — > D*^ l~ v\ These model dependent extractions have unknown theoretical errors. 
Following the model dependent extractions we describe the measurement of |V^b|jF(l) from 
the differential decay distribution dI^{B ^ D* i u)/dy. Using a model prediction for the 
corrections to the heavy quark symmetry limit at y = 1 we then present a value for \Vcb\- 
We also compare the results for the shape of J^{y) with model predictions. 



A. Model dependent determinations 

There are several quark models that predict the normalization and dependence of the 
form factors describing the decay B ^ D* i u. When extracting model dependent values 
for \Vcb\ based on predictions of the form factors and their dependence it is important 
to determine the efficiency for the model under investigation because the acceptance can 
vary with different values of Ai, A2/A1 and V/Ai. The form factor ratios and their 
dependence for the different models considered here are in Table 0. Using the partial width 
from Eq. (PDj), \Vcb\ can be calculated from model predictions of the rate for B ^ D* £ u. 
The calculation includes a correction for the difference in the efficiency predicted with each 
model; the values for \Vcb\ obtained with four different models are given in Table |V11| . 



B. \Vcb\ determined from the dr{B —>■ D* £ y)/dy distribution 

In this section the extraction of from the dT{B —>■ D* [ i')/dy distribution is 

described. Three fits are performed. Using the assumption that /+_ = /oo = 0.5, we first 
extract |V"c6|J?^(l) separately from the differential decay distributions for B^ D*^ £~ v 
and for B~ D*^ £~ v events. We then simultaneously fit the B^ — ^ D*^ £'^ v and the 
B~ — >■ D*^ £~ v events, which results in a determination of |V^f,|jF(l) that is independent of 
/oo/ /+-• 

We use the unbinned maximum likelihood method developed for fitting the multidi- 
mensional decay distribution oi D ^ K*W events I^Tj. Application to the extraction of 
\Vcb\J^{y) is straightforward since only one dimension is used in the fit. It is convenient 
to make the likelihood a function of Vcb by normalizing our probability distribution func- 
tion to the number of observed events. This requires the inclusion of a Poisson probability 
factor 



D* and lepton pairs are selected in the Signal Region of the C vs. MM"^ plane as 
already described. However, since we are now performing an unbinned likelihood fit, we 
apply a M^n cut as well as a 6^ cut. The M^-k cuts are \Mktt — mDo\ < 25 MeV and 
iMxn — mDo\ < 20 MeV for B^ D*^ £~ v and B^ D*^ £~ v respectively, where m^jo 
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is the nominal mass. In order to have low background levels in the unbinned likelihood 
fit, the tighter B' D*^ t v cut was chosen to optimize signal to background as much as 
possible, without introducing a systematic bias from reproduction of signal shapes by the 
Monte Carlo simulation. 

Eq. ( p^ ) can be rewritten as 

^ = ^(2/)|K5r-^'(y), (33) 

where all the known quantities have been folded into a single function Q[y)- jVcbl and T[y) 
are to be determined by fitting the data. Since we want to determine the product \ych\T{y) 
at y = 1, we approximate the unknown function T{y) with an expansion about ?/ = 1, 

T{y) = ^(1) [l - a'{y - 1) + b{y - if] , (34) 

thus defining the fit parameters to be |V^b|JF(l), and b. We use rather than a because 
the B and D* meson wave functions have maximum overlap at ?/ = 1, so the first derivative 
of J-'{y) must be negative at that point. Our results for and b can be used to compare 
the shape of J-'{y) determined from data with theoretical model predictions ||53| . 

The variable y is the D* energy to mass ratio in the rest frame of the decaying B. 
However, in data we measure E£)*/m£)* in the laboratory frame, symbolized by y. This 
quantity is a function of y and the B meson momentum p^, but since the direction of 
is unknown, the dN/dy distribution of the data is not directly observed. We therefore fit 
dN/dy, and naturally incorporate the smearing due to the B motion in the same way as 
the smearing due to detector resolution (which is much smaller). This is accomplished by 
a function 7l{y,y), which is convoluted with G{y)J-'{y) in order to arrive at a probability 
distribution function that can be compared with the data, 

ryo 

G{y,a,b)=J^ giy)J^\y)n{y,y)dy, (35) 

where yo is the upper kinematic limit of y, given by = 0. The resolution function TZ{y, y) is 
determined from Monte Carlo simulation, and allows G{y, a, b) to correctly include changes in 
detection efficiency with variations in the parameters of JF, thus reducing model dependence. 
The full likelihood function is 

£(|K6|,a,&) = e-^(l^-l''^''')-^''n fiV(|K,|,a,6)G(y,,a,6) , (36) 

1=1 \ b J 

with 

/yo 
e{y)G{y):F\y)dy (37) 

Nb = ricomb + ?T-D"«P + ^uncorr (38) 

^uncorr -Puncorr (l/i) • (39) 

b 

N{\Vcb\,a,b) is the expected number of reconstructed signal events, where eiy) is the effi- 
ciency for detecting B ^ D* £ u events determined from the Monte Carlo, tb is the B 
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meson lifetime and / (meaning either or /oo), iVT(4S), ^d* and Bdo have been defined 
in Eqs. § and ^. The total number of background events passing our cuts is Nf,. The 
exponential factor in Eq. (|36| ) comes from the Poisson probability mentioned earlier. The 
backgrounds are represented by Eqs. (^) and (|39D . They have been divided into three 
types, combinatoric, B D* X i u, and uncorrelated background, and their normalization 
and y dependence is represented by Eq. (pOl). The combinatoric background yield and y 
distribution, (ncomb-Pcomb(^)), are obtained from sidebands in the 6m distribution. The cor- 
related background yield, nD**iu, is obtained from Table | or H and its y distribution, 
PD**w{yi)i is taken from Monte Carlo generated as in Section |V1 C| according to the ISGW 
model. The uncorrelated background yield and y distribution, nuncorr-Puncorr(y) is obtained 
from the measured inclusive B D* and B ^ spectra, assuming an isotropic angular 
distribution between the uncorrelated leptons and D*^s. The three y distributions of the 
backgrounds to 5° — D*~^ i~ v and B~ D*^ i~ v are shown in Fig. 

The data yields as a function of y for B^ D*^ £^ v and i?" D*^ i~ v events 
are shown in Fig. [1^. The values for |V^fe|jF(l) and from linear fits to these data are 
given in Table |V111| . We have fixed h to zero in these fits because the statistics of the 
individual modes are not sufficient to constrain the curvature of The first error listed 

in Table |Vlli| is statistical and the second error is an estimate of the systematic error, 
including the uncertainty in the B lifetimes. The statistical uncertainty in the backgrounds 
has been taken into account by varying the background levels in the fit by ±1 standard 
deviation of their measured central values, and is included in the statistical error. 

The contributions to the systematic error on | V^6|jF(l) and are listed in Table |X|. The 
errors on the absolute efficiencies are the same as for the branching fraction measurements. 
The errors that are detailed in Table IX are specific to the dV/dy fits and account for the 
uncertainty in the efficiency as a function of y and uncertainties in the fitting procedure. As 
a check of the fitting technique we have performed a binned fit to the efficiency corrected 
dr/dy distribution, using an analytically smeared form 



27P ^Jy^^n{y) VI/ - 1 

where 7/? is the known momentum to mass ratio of the decaying B meson and (ymin,ymax) is 
the range of possible values of y given a value of y (see Fig. H). The results agree well with 
the results obtained from the unbinned likelihood fit. The systematic uncertainties due to 
the fitting procedure were estimated by comparing the results of the unbinned and binned 
methods in 36 Monte Carlo samples, each with the same statistics as our data sample. 

We have also fit the differential decay distribution of the two decays modes B^ — > 
D*'^ i~ V and B~ —>■ D*^ i~ v simultaneously by calculating the likelihood of the two 
modes separately and maximizing their product. The result of the simultaneous fit is shown 
in Fig. |TT], along with the combined data and background level for the B^ —>■ D*^ i~ v 
and B^ D*^ i~ v modes. From a linear fit to these data we obtain the most precise 
measurement to date of |V^b|jF(l) and [p^, 

\V^^\T{\) = 0.0351 ± 0.0019 ± 0.0018 ± 0.0008 (41) 

= 0.84 ±0.12 ±0.08, (42) 

independent of /+_//oo; where the first error is statistical, the second is systematic, and the 



third is due to the input lifetimes [56|. With the statistics of both modes combined we can 
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also perform a three parameter fit, wliere tlie quadratic term in tlie expansion for J^{y) is 
retained. The results of this fit are shown in the last row of Table [VIII| . For comparison, we 
provide in Table ^ the a? (and h) values predicted by the different decay models. 

To extract a value of \Vch\ one requires a prediction for ^^-"(1). Using J^il) = 0.97 ± 



0.04 11,13 we obtain, 

\Vcb\ = 0.0362 ± 0.0019 ± 0.0020 ± 0.0014, (43) 
where the third error is the quoted theoretical uncertainty in J^(l) and the uncertainty 



in the B lifetimes has been included in the systematic (second) error [^. This result is 
consistent with previously published values of \Vcb\ (after scaling for differences in D and D* 
branching fractions and B lifetimes) from both inclusive and exclusive analyses, although 
the unknown model dependence of the previous \Vcb\ extractions makes detailed comparisons 
difficult. Model dependent values of \Vcb\ in Table |V11| are also consistent with this result. 

Other recent estimates are J^(l) = 0.96 ± 0.03 fT^ and J^(l) < 0.94 The value 

of J^{1) used to obtain Eq. (^) is consistent with both of these results. However, Ref. 
Il58| also claims an "educated guess" of J^(l) = 0.89 ± 0.03, which would lead to \Vcb\ = 
0.0395 ± 0.0021 ± 0.0022 ± 0.0012. 

The product \Vcb\J^{y) as determined from our fits to the combined B^ — > D*^ £^ v and 
5" D*^ i~ V data is shown in Fig. [l^. Fig. 0(a) shows the result obtained by keeping 
only the linear term in the expansion of T[y\ while (b) shows the result for a second order 
polynomial. The dotted lines show the contours for la variations of the fit parameters. All 
the lines are functions of the variable ?/, which is equal to the D* energy to mass ratio, 
E]j*/mj:,*, evaluated in the B rest frame. The points are the data. As discussed before, 
we do not measure y but rather y, which is Ej:i*/m£,* in the laboratory frame. We have 
therefore binned the data in the average of the maximum and minimum possible values 
of y for each measured y. This average is denoted by yAijj)- Unfortunately yAijj) is not 
an unbiased estimator or ?/, and although yAijj) has the same kinematic limits as the 
smearing between these two variables is significant compared to the bin size. Therefore, the 
data shown this way do not exactly correspond to the product \Vcb\^{y)- Note, however, 
that the lines shown are not the result of a fit to these data points, but were obtained 
from our unbinned fit which correctly accounts for the boost between the B frame and the 
laboratory frame. 

The shape of J^{y) determined from the fit to the data can be used to obtain the functional 
form of dr{B D* i u)/dq'^] (where is the true of the decay, not the smeared variable 
measurable experimentally) which is useful to test hypotheses of factorization in B meson 
decay . Values of dr{B D* i v)/ dq^ as determined from the data are given in Table ^ 



for commonly used q points. 



VIII. CONCLUSIONS 

Using the CLEO II data sample we have reconstructed both B^ D*~^ i~ v and B~ — >■ 
V decay modes in the D* decay chains: D*+ D^7i~^ and D*^ D^n° with 
— > K^TT^. We have measured the exclusive branching fractions. 
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B{B^ D*+ r u) = (0.5//oo)[4.49 ± 0.32(stat.) ± 0.39(sys.)]% (44) 
B{B- D*^ r v) = (0.5//+_)[5.13±0.54(stat.) ±0.64(sys.)]%, (45) 

which depend on the relative production of charged and neutral B mesons at the T(4S') 
resonance. We have also presented these results in a way that requires no knowledge of 
Z+./Zoo, but does depend on B lifetime measurements from other experiments. With the 
assumption of equal partial widths, r{B^ D*~^ £~ v) = r{B~ D*^ i~ u), we have 
determined the product 

= 1.14 lb 0.14(stat.) ± 0.13(sys.), (46) 

Joo tbo 



and by using existing B lifetime ratio measurements we have solved for /+-//00 = 
1.04 ± 0.13(stat.) ± 0.12(sys.) ± 0.10 (lifetime). This measurement verifies at the 10% level 
the assumption that /+_ = /qo, widely used to calculate B meson branching fractions with 
data collected at the T(4S') resonance. Also with the assumptions of equal partial widths 
and /+_ + /oo = 1, and using external B lifetime measurements, we have determined the 
partial width 

T{B ^ D* iu) = [29.9 ± 1.9(stat.) ± 2.7(sys.) ± 2.0(lifetime)]ns-\ (47) 

independent of Z+./Zoo. From this width we have determined the B ^ D* £ u branching 
fractions with dependence only on the charged to neutral B lifetime ratio (but not the 
individual hfetimes), 

13{B- r u) = 1.1 X i3(5° ^ D*+ r u) 

= [5.03 ± 0.32(stat.) ± 0.45(sys.) ± 0.24(lifetime)]%. (48) 

where the factor 1.1 carries no error because we are assuming that the ratio of branching 
fractions exactly equals the value used for the lifetime ratio (the uncertainty in our knowledge 
of the lifetime ratio results in the third error). 

Taking advantage of theoretical constraints on the normalization and dependences 
of the form factors for the B D* i u decay provided by HQET, a combined fit to 
the differential decay distributions for B^ D*~^ i~ v and B~ — >■ D*^ i~ v results in a 
determination of |K:{,|JF(1), 

\Vd\T{\) = 0.0351 ± 0.0019(stat.) ± 0.0018(sys.) ± 0.0008(lifetime), (49) 

which is also independent of /+_//oo- The explicit dependence on the B^ and B^ lifetimes 
is given in Ref. [Q. The form factor normalization J-'{1) is believed to be predictable 
with small theoretical uncertainty, and this permits a precise determination of the CKM 
matrix element. Using the prediction of Ref. for the normalization of J-', we obtain 

\Vcb\ = 0.0362 ± 0.0019(stat.) ± 0.0020(sys.) ± 0.0014(model), where the uncertainty in the 
B lifetimes is included in the systematic error. 
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TABLES 



TABLE L A summary of the events in the B D*^ v Signal Region and Correlated 
Background Region. In the top 2 rows the first error is statistical and the second is the systematic 
uncertainty due to the choice of functions used to fit the data. The bottom two rows have been 
computed using Eq. ( |To| ) and the uncertainties in and Tc have been propagated into the systematic 
error. The numbers in parenthesis are not independent measurements, but merely functions of the 
other numbers in the bottom two rows. 



D*+l- 


Signal 
Region 


Correlated 
Background Region 


5m Signal 
Scaled 5m Sideband 
Uncorrelated Background 

Continuum 
Correlated Lepton Fakes 
B D*Tu and B D*Ds 


457 ± 23 ± 9 
47 ± 6 ± 1 
5.0 ±0.5 
5±7 
2± 1 


70 ± 11 ± 1 
9±5 
19 ± 2 
0±6 
1 ± 1 
4±3 


Net Yield [Ns and Nc) 


398 ± 25 ± 9 


37± 14± 1 


Ns{B° D*+ £- u) 
Nc{B^ D*+ £- u) 


376 ± 27 ± 16 
(22 ±9 ±11) 


(6 ±0.5) 
31 ± 14± 1 


TABLE II. A summary of the events in the B~ D*^ v Signal Region and Correlated 
Background Region. In the top 2 rows the first error is statistical and the second is the systematic 
uncertainty due to the choice of functions used to fit the data. The bottom two rows have been 
computed using Eq. (|l^) and the uncertainties in and have been propagated into the systematic 
error. The numbers in parenthesis are not independent measurements, but merely functions of the 
other numbers in the bottom two rows. 




Signal 
Region 


Correlated 
Background Region 


5m Signal 
Scaled 5m Sideband 
Uncorrelated Background 

Continuum 
Correlated Lepton Fakes 
B D*Tv and B D*Ds 


476 ± 25 ± 10 
144 ± 10 ± 3 

8±1 

6±7 

2±2 


94± 15± 1 
34 ±7 
23 ±3 
9 ±11 
1 ± 1 
2±2 


Net Yield {Ns and Nc) 


316 ± 28 ± 10 


25 ± 20 ± 1 


N,{B- D*^ £- v) 
Nc{B- L>*° r 9) 


302 ± 32 ± 13 
(14 ± 14 ± 7) 


(5 ±0.5) 
20 ± 20 ± 1 
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TABLE III. Contributions to the fractional errors of the detection efficiency for the 
Z)*+ V and B~ — > D*^ l~ v decays. Errors common to both modes are entered only once 
in the column labeled "Both" . 



hource 


Aes/es (%) 




D — > L) ^ t V 


rjotn 


D — > JJ t V 


Form Factors 


3.0 




byji Sideband Normalization 0.8 




1.7 


Variation of Cuts 2.3 




1.0 


Lepton Efficiency 


2.2 




K~t:+ Efficiency 


4.0 * 




Slow 7r+ Efficiency 


5.0 




Slow -K^ /-K^ Efficiency Ratio - 




7.0 


Total Uncorrelated 2.3 




7.3 


Total Correlated 


7.4 




* This error is correlated with the systematic error of the CLEO II 


i, ZV^ iT~^ Ttt Q n oil m (T 

^ 1\ /I Ui cLliL/lilli^ 


fraction. 






TABLE IV. Results for product branching fractions. 




Product Branching Fraction 




Result xlO^ 


e(T(45) ^ B^B^)B{B^ D*+ r i?)B{D*+ D^n 


^)B{D^ ^ K-7T+) 


6.0 ±0.43 ±0.55 


13{T{4S) B+B-)B{B^ D*+ £- v)B{D*^ D^tt 


^)B{D^ ^ K-^+) 


6.4 ±0.68 ±0.73 


B{B D*+ r u)B{D*+ D^^+)B{D^ K-7r+) 




1.6 ±0.7 ±0.3 


B{B D*o X £- i?)B{D*° D^tt^)B{D° K-^+) 




1.4 ± 1.5 ±0.3 



TABLE V. Comparison with previously published results for B{B^ D*^ £^ v) and 
B{B~ D*^ i~ v). All previous results have been rescaled to use the CLEO II and D* 
branching ratios (except for the CLE01.5 B~ number which depends nonlinearly on these branch- 
ing fractions), and use /-| = /oo = 0.5. 

Experiment B{B^ D*+ £- u) (%) B{B- r i?) (%) 

CLEO II {tb+/tbo = 1.10 ± 0.11) 4.57 ±0.29 ±0.41 ±0.22 * 5.03 ± 0.32 ± 0.45 ± 0.24 * 
CLEO II (/+_ = /oo = 0.5) 4.49 ±0.32 ±0.39 5.13 ± 0.54 ± 0.64 

ARGUS gl 4.9 ±0.5 ±0.6 6.4 ±1.5 ±1.4 

ARGUS (g 4.5 ± 0.3 ± 0.4 

CLEO 1.5 IJ3 4.1 ±0.5 ±0.6 4.1 ±0.8 ±0.9 

* When tb+/tbo = 1.10 is used to extract the CLEO II branching fractions, all errors for 
the two different exclusive modes are fully correlated. 
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TABLE VI. Predictions for the form factor ratios and their or y dependence, and the predic- 
tions and dependence for the common form factor A\. The pole forms have been abbreviated as: 
Pi = 1 - g2/6.342 and P2 = 1 - gV6-732 where q is in GeV. £{(^) stands for exp[-0.03(g^„^ - q^)\. 
For the form factor ratios of the Neubert model [^0| A^(y) = 2.b/{y + 1). 



ISGW H BSW g KS g 



Neubert Ol 



A2{q^)/A^{q^) 
Ai{q^) 



1.27 
1.14 

QM8{q^) 



imP2/Pi 

1.06 

0.65/P2 



1.00/Pi Af(y)[1.35 - 0.22(y - 1) + 0.09(y - 1)^] 
1.00/Pi Af(y)[0.79 + 0.15(y - 1) - 0.04(y - 1)^] 
0.70/Pi 0.86[2/(y + l)]^-^ 



TABLE VIL Model dependent predictions of the B ^ D* I v partial width, detection efficien- 
cies for each model, and \ Vcb\ values derived from the measured partial width. An additional 3.5% 
systematic error due to B lifetime measurements is common to all the \Vcb\ values given. 



ISGW 


31 


1 BSW 


41 


1 KS 


42] 


Neubert 


39| 


T{B ^ D* I v)\V,b\'^ 


24.6 ps 


-1 




21.9 ps-i 


25.8 ps-i 


29.0 ps" 


1 



emodel/eNeubert 1-00 0.97 0.98 1 

iKbl X 10^ 34.8 ±1.1 ± 1.6 37.5 ±1.2 ± 1.7 34.4 ±1.1 ± 1.5 32.2 ±1.0 ±1.4 



TABLE VIII. Values for |Vcfe|.F(l), a? and b determined by the fits to the dV/dy distributions. 
The first error is statistical and the second systematic (including B lifetimes). The values for the 

individual modes, B° D*+ i' 9 and B' D*^ i' u, assume /oo = f-\ = 0.5. The results of 

the combined fit (bottom two rows) are independent of this assumption. 



Mode 


\v,b\H^) X io3 




b 




34.7 ± 2.5 ± 1.8 


0.80 ±0.17 ±0.08 


0.0 


B- P>*o £- u 


35.7 ±2.8 ±2.4 


0.94 ±0.20 ±0.08 


0.0 


B ^ D* e 9 


35.1 ± 1.9 ± 1.9 


0.84 ±0.13 ±0.08 


0.0 




35.3 ± 3.2 ± 3.0 


0.92 ±0.64 ±0.40 


0.15 ± 1.24 ±0.90 
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TABLE IX. Estimates of the systematic error on Vcb and determined from the fits to the 
dT/dy distribution. The fractional error on Vcb is the same as the fractional error on the product 
I V^ft|.7-"(1) (the uncertainty in theoretical estimates of .F(l) is not included in this table). The 
error on the absolute efficiency is taken from Table III, excluding the 4% contribution from Kit 
reconstruction, which is also excluded from the error for B£)0. 



source 


^0 ^ 


- £' u 


B- D*^ 




B 


D* £u 








AVcb/Vcb 




^Vcb/Vcb 




B lifetime 


2.9% 




3.6% 




2.3% 




Absolute Efficiency 


3.3% 




4.8% 




3.3% 






1.9% 




3.4% 




2.1% 




Slow TT efficiency shape 


1.2% 


3.5% 






0.8% 


2.4% 


Fitting systematics 


3% 


9% 


3% 


9% 


3% 


9% 


Total 


5.8% 


9.7% 


7.5% 


9.0% 


5.5% 


9.3% 



TABLE X. Model predictions for the form factor parameters a} and h. These values were 
determined by fitting the dL / dy distributions predicted by each model. Results are given with and 
without h fixed to zero. 



Parameter 


This Exp. 


ISGW 131] 


BSW |41| 


KS @ 


Neubert ^ 


a2 (h = 0) 


0.84 ±0.13 ±0.08 


0.91 


0.77 


0.83 


0.48 




0.92 ± 0.64 ± 0.40 


0.88 


0.84 


1.10 


0.59 


b 


0.15 ±1.24 ±0.90 


-0.1 


0.1 


0.6 


0.28 



TABLE XL Values for the differential decay rate at selected points. These were calculated 
from Eqs. ( |32|) and (||) with the results of the maximum likelihood fit given in Eqs. @) and (||). 
The errors are combined statistical and systematic. 





2 


2 






dV/dq^ (ns-iGeV-2) 


1.64 ±0.23 


1.90 ±0.25 


2.32 ±0.31 


3.17 ±0.44 
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FIGURES 




FIG. 1. Kinematic boundary for B — ^ D* £ v decays in the plane of C vs. MM"^, for lepton 
momentum in the range 1.4 < |p^| < 2.4 GeV (sohd line). This is the Signal Region for this 

analysis. The dots, including those outside the kinematic boundary, are Monte Carlo signal events. 
Final state radiation and bremsstrahlung occasionally force reconstructed events outside of the 
kinematic boundary. 
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FIG. 2. Mass difference {5m) distributions of events falling in the Signal Region of the decay 
modes (a) D*^ t~ v and (b) B~ D*^ £~ v. All candidates are required to have Mk-k 

within 100 MeV of the measured mass. 
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FIG. 3. Invariant mass of K vr"*" combinations passing a 5m cut in the Signal Region of the 
decay modes (a) S° D*+ fr v and (b) D*^ £- v. The dashed histo grams show the scaled 

Mk-k distributions of events in the Signal Region, but in the 5m sideband. The solid lines show 
the fits used to determine the yields. 
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FIG. 4. Distribution of Monte Carlo B D** i u events in the plane of C vs. MM^, (a) for 
lepton momentum in the range 1.4 < |p^| < 2.4 GcV and (b) 0.8 < \pi\ < 1.4 GcV. The kinematic 
boundary for B D* I v decays is also shown in each plot (solid line). The dashed line in the 
right plot indicates the boundary of the Correlated Background Region. 
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FIG. 5. Invariant mass of K vr"'" combinations passing a 6^ cut in tli6 Correlated Background 
Region of the decay modes (a) v and (b) — > D*^ l~ v. The dashed histograms 

show the scaled Mk-k distributions of events in the bm sideband. The solid lines show the fits used 
to determine the yields. 
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FIG. 6. Monte Carlo reconstruction efficiency versus momentum (filled diamonds) for charged 
slow pion candidates. Overlaid (solid curve) is the efficiency as determined from the inclusive D*+ 
decay angle distribution in data. The la variations in the parameters of this curve are shown by 
the dotted lines. The dashed line indicates the minimum momentum cut used in this analysis. The 
efficiency shown here includes the geometric acceptance of 71%. 
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FIG. 7. Monte Carlo reconstruction efficiency versus momentum (filled diamonds) for neutral 
slow pion candidates. Overlaid (crosses) are the results of the r/ and studies which compare 
the ratio of neutral to charged pion efficiency in data and Monte Carlo. These two points were 
calculated assuming that the Monte Carlo exactly simulates the charged pion efficiency. The 
efficiency shown here includes the geometric acceptance (not constant with 7r° momentum) , which 
accounts for much of the inefficiency. 
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FIG. 8. The smearing between the variables y and y for a lepton in the momentum range 1.4 



to 2.4 GeV. The points are from Monte Carlo events generated according to the model of Ref. |39|. 
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(a) Backgrounds for B°^D * 1" v 
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FIG. 9. The y distributions of the backgrounds for (a) D*^ fr v events and for (b) 

B~ D*^ (r u events. The sohd hne is the combinatoric background, the dashed line is the 
correlated background, and the dotted line is the uncorrelated background. The area of the curves 
are normalized to represent the background levels in the data. 
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FIG. 10. Differential yield, clN/dy, for (a) 5° D*+ v events and (b) B' L>*° v 
events in the data, with the projections of the unbinned fits superimposed. The solid histogram 
represents the result of the fit using J^{y) = 1 — a^(j/ — 1). The dashed histogram shows the level 
of the background from all sources. 
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FIG. 11. Differential yield, dN/dy, for B' ^ D 
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combined with the projection of the unbinned simultaneous fit superimposed. The solid histogram 



represents the result of the fit using J^{y) 
of the background levels. 



1 — a^(y — 1). The dashed histogram shows the sum 
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FIG. 12. The product |V^fe|^(y) (solid line) as determined by our fits fits to the combined 
D*^l^ and 0*^1" data, using (a) a linear expansion and (b) a quadratic expansion of the form 
factor J-{y). The dotted lines show the contours for la variations of the fit parameters. The 
points are data for the square root of the measured decay rate divided by the factors other than 
J-{y) and Vcb of Eq. (^) (error bars are statistical only). As explained in the text, the data are 
binned in yA{y), which is not an unbiased estimator of y. Therefore, the data points do not exactly 
correspond to the product \Vcb\J^{y)- 
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